Early Dark Energy at High Redshifts: Status and Perspectives 
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Early dark energy models, for which the contribution to the dark energy density at high redshifts 
is not negligible, influence the growth of cosmic structures and could leave observable signatures 
that are different from the standard cosmological constant cold dark matter (ACDM) model. In 
this paper, we present updated constraints on early dark energy using geometrical and dynamical 
probes. From WMAP five-year data, baryon acoustic oscillations and type la supernovae luminosity 
distances, we obtain an upper limit of the dark energy density at the last scattering surface (Iss), 
f2EDE(ziss) < 2.3 X 10~^ (95% C.L.). When we include higher redshift observational probes, such as 
measurements of the linear growth factors. Gamma- Ray Bursts (GRBs) and Lyman-a forest (Lya), 
this limit improves significantly and becomes nEDB(-Ziss) < 1.4 x 10~^ (95% C.L.). Furthermore, we 
find that future measurements, based on the Alcock-Paczyhski test using the 21cm neutral hydrogen 
line, on GRBs and on the Lya forest, could constrain the behavior of the dark energy component 
and distinguish at a high confidence level between early dark energy models and pure ACDM. In this 
case, the constraints on the amount of early dark energy at the last scattering surface improve by a 
factor ten, when compared to present constraints. We also discuss the impact on the parameter 7, 
the growth rate index, which describes the growth of structures in standard and in modified gravity 
models. 



I. INTRODUCTION 

Current cosmological observations, such as the cosmic 
microwave background (CMB) measurements of temper- 
ature anisotropies and polarization [l| and the redshift- 
distance measurements of Type la Supernovae (SNIa) at 
z < 2 0], have demonstrated that the Universe is now 
undergoing an accelerated phase of expansion and that 
its total energy budget is dominated by the dark energy 
component. The nature of dark energy is one of the 
biggest unsolved problems in modern physics and has 
been extensively investigated in recent years, both un- 
der the theoretical and the observational point of view. 
The simplest candidate of dark energy is the cosmologi- 
cal constant (ACDM): however, this model suffers of the 
fine-tuning and coincidence problems 0]. To lift these 
tensions, many alternative dynamical dark energy mod- 
els, such as quintessence [3l, phantom K-essence @ 
and quintom Q , have been proposed (for a review see 
Ref.d). 

Among all the dynamical dark energy models, we will 
focus here on early dark energy (EDE) ones, in which 
a small fraction of dark energy is present up to the 
last scattering surface (Iss), unlike ACDM for which 
^^de(21ss) — 0, The differences between early dark en- 
ergy models and pure ACDM are particularly evident at 
high redshifts, over a large fraction of the cosmic time, 
when the first structures form. EDE has been shown to 
influence the growth of cosmic structures (both in the 
linear and in the non-linear regime), to change the age 



of the universe, to have an influence on CMB physics, 
to impact on the reionization history of the universe, to 
modify the statistics of gian t arcs in str ong cluster Icnsing 

statistics I, M, M MM, M M M m- 

The reason for addressing this particular model of dark 
energy is also driven by the increasing availability of high 
redshift observations that somewhat bridge the gap be- 
tween the CMB and very local cosmological probes. Dark 
energy studies rely mainly on: the clustering properties 
of luminous red galaxies [l^, baryonic acoustic oscilla- 
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tions (BAG) 
redshift SNIa 

as standard candles, and also the number of high redshift 
galaxies and clusters of galaxies [1^ that could be studied 
with deep field observations. Of course, since the present 
data sets seem to be in good agreement with pure ACDM 
model, we hope that these observations could open up a 
new high redshift window on the properties of dark en- 
ergy and possible allow to confirm or disproof its redshift 
evolution. Such a difficult measurement would have pro- 
found implications for physical cosmology and particle 
physics. 

In this paper we extend our studies of dark energy 
behavior into the redshift range 2 < z < 1100, which 
covers the so-called dark ages, and present constraints 
on EDE from current and future cosmological observa- 
tions. We will use CMB, BAG, SNIa, Lya forest ob- 
servations of growth factors and matter power spectrum 
and the GRBs from current available data sets, while for 
forecasts we will add to Planck-like observations, some 
higher redshift GRBs, some improved Lya constraints 
and a measurement that could be potentially very im- 
portant of the so-called Alcock-Paczyhski (AP, [2^) test 
using 21cm maps [13, HH- As we will see, the partic- 
ular parameterization of EDE chosen will be such that 
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the higher the redshift the tighter the constraints will be 
on the parameters describing the given model. Thereby, 
futuristic intergalactic medium Lya forest data, GRBs 
distance moduli and the 21cm maps, before or around 
reionization, are expected to be very promising probes of 
dynamical dark energy models (see also Ref. [l3| for the 
impact of these models on the reionization history of the 
universe). We note that the behavior of dark energy in 
this redshift range has been exploited by Ref.js^, who 
constrained the EDE density of to be less than 2% of the 
total energy density. Here, we will improve this limit by 
adding several different observations: as mentioned be- 
fore this is not only important to get stronger constraints 
on some cosmological parameters but also to understand 
the consistency of different probes and how systematic 
effects irnpact on the final derived measurements (e.g. 
Rcfs.jH, [131). Of course, we stress that while our ap- 
proach has the great advantage of using dynamical and 
geometrical cosmological probes in the linear or quasi- 
linear regime, other approaches, based for example on 
haloes concentrations, could also be envisaged and must 
rely on an accurate comparison with numerical simula- 
tions [l6j . 

Our paper is organized as follows: in Sec. II we de- 
scribe the theoretical framework of the early dark energy 
model and the datascts we used. Sec. Ill contains our 
main global fitting results from the current observations. 
In Sec. IV we present the forecasts for the future mea- 
surements while Sec. V is dedicated to the conclusions 
and discussion. 



II. METHOD AND DATA 
A. Parameterization of early dark energy 

We decided to use the mocker model introduced in 
Ref. [131 motivated by the following two observational 
facts: i) the best-fit model is pure ACDM model and 
the amount of dark energy at the last scattering sur- 
face is constrained to be close to zero from CMB and big 
bang nucleosynthesis (BBN) observations; ii) an equation 
of state of the dark energy component which is rapidly 
evolving (i.e. dw/dz 3> 1) seems to be ruled out by high- 
redshift SNIa [H. 

For these EDE models the parameterization reads: 



li'EDE(a) = -1 



1 - 



1 



(1) 



where a = 1/(1 -I- z) is the scale factor, wq is the present 
equation-of-state of dark energy and C characterizes the 
"running" of the equation of state. Consequently, the 
evolution of dark energy density can easily be obtained 
via energy conservation as: 
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Figure 1: Left panel; the evolution of dark energy density 
and equation of state for three models, too = —0.972 and C = 
1.858 (EDEl), Wo = -0.95 and C = 2.5 (EDE2) and pure 
ACDM. Right panel: the linear growth factor G{a) = S{a)/a 
normalized at 2 = oo for the three models. The shaded area 
represents the redshifts that we will be mainly probing in the 
present work (2^2 — 20). 



In Fig. [T] we plot the dark energy density (upper part 
of the left panel) and equation of state (bottom part of 
the left panel) as a function of redshift for three differ- 
ent models: pure ACDM (dotted blue hue), EDEl (solid 
black line) which has iwo,C) (-0.972,1.858) and 
EDE2 (dashed red line) which has (wq, C) = (-0.95, 2.5). 
All these models fit the CMB and the lower redshift SNIa 
constraints very well. In the right panel, we also show the 
linear growth rate (growth factors divided by the scale 
factor) normalized at z = cx) and compare the different 
evolutions with redshifts in the three cases: differences of 
the order of 20% are visible at z ~ 3. In the two panels 
we also show as a shaded vertical band the region in the 
redshift range z = 2 — 20, a period which is roughly 25% 
of the age of the universe. This is the redshift range we 
will be (mainly) focussing on in the rest of the paper in 
order to discriminate between the different dark energy 
models. 

We stress that this is just one of the possible paranie- 
terizations for early dark energy models, another one has 
been suggested by Ref. 35 1 and recently used in Ref. [l7| . 
However, we prefer to use the mocker model in order to 
compare with the findings of Ref.^^ and because this 
parameterization has a smooth redshift derivative at low 
z for w{z). We note that one of the most important pa- 
rameters is the amount of dark energy during the struc- 



ture formation period and this is given by {a^q is the 
matter-radiation equality scale factor): 



(2) 
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and we will also quote this value in the rest of the paper 
in order to compare with other works and constraints as 
well (e.g. Refs.[i,[il). 

We will also use some growth factors measurements 
obtained mainly through Lya forest observations and 
one lower redshift from Sloan Digital Sky Survey (SDSS) 
galaxies. For this purpose we will compute the growth 
factors for the EDE models with a subroutine imple- 
mented in the publicly available Monte Carlo Markov 
Chains (MCMC) package CosmoMC ^ that solves the 
second order differential equations for the growth factor 
as in Ref.[l3|. Furthermore, we also modify CosmoMC 
to include the perturbations of dynamical dark energy 
models generally as done in Refs.|37|. 



B. Current and Future Datasets 

1. Current Observations 

We will rely here on the following cosmological probes: 
i) CMB anisotropics and polarization; ii) baryonic acous- 
tic oscillations in the galaxy power spectra; iii) SNIa dis- 
tance moduli; iv) GRBs distance moduli; v) the Lya for- 
est growth factors and matter power spectrum measure- 
ments. 

In the computation of CMB power spectra we have 
included the WMAP five-year (WMAP5) temperature 
and polarization power spectra with the routines for 
computing the likelihood supplied by the WMAP team 
jssl l . Besides the WMAP5 information, we also use some 
distance-scale indicators. 

BAOs have been detected in the current galaxy redshift 
survey data from the SDSS and the Two-degree Field 
Galaxy Redshift Survey (2dFGRS) [H, HO, Hi 113 • The 
BAO can directly measure not only the angular diameter 
distance, Dj^(z), but also the expansion rate of the uni- 
verse, H{z), which is powerful for studying dark energy 
[111 . Since current BAO data are not accurate enough for 
extracting the information oi Da{z) and H{z) separately 
[4^ , one can only determine an effective distance [T^ : 



{l + zfD\{z) 



H{z) 



1/3 



(4) 



In this paper we use the Gaussian priors on the distance 
ratios rs{zd)/Dy{z): 



rs{zd)/D,{z = 0.20) 
r,(zd)/D,(z = 0.35) 



0.1980 ±0.0058 , 
0.1094 ± 0.0033 , (5) 



with a correlation coefficient of 0.39, extracted from the 
SDSS and 2dFGRS surveys [i^l, where Vg is the comoving 
sound horizon size and Zd is drag epoch at which baryons 
were released from photons given by Ref. [i^ . 

The SNIa data provide the luminosity distance as a 
function of redshift which is also a very powerful mea- 
surement of dark energy evolution. The supernovae data 



we use in this paper arc the recently released Union com- 
pilation (307 samples) from the Supernova Cosmology 
project which include the recent samples of SNIa from 
the (Supernovae Legacy Survey) SNLS and ESSENCE 
survey, as well as some older data sets, and span the red- 
shift range < z < 1.55. In the calculation of the likeli- 
hood from SNIa we have marginalized over the nuisance 
parameter as done in [4^ and ignored the systematic er- 
rors to improve our results^. 

Furthermore, we make use of the Hubble Space Tele- 
scope (HST) measurement of the Hubble parameter 
H(j = 100 /i km s^^ Mpc^^ by a Gaussian likelihood 
function centered around h = 0.72 and with a standard 
deviation cr = 0.08 [H. 

In order to constrain the early dark energy model at 
high redshifts, we also include the GRBs, linear growth 
factors and Lya forest data. 

GRBs can potentially be used to measure the luminos- 
ity distance out to higher redshift than SNIa. Recently, 
several empirical correlations between GRB observables 
were reported, and these findings have triggered inten- 
sive studies on the possibility of using GRBs as cosmo- 
logical "standard" candles. However, due to the lack of 
low-redshift long GRBs data to calibrate these relations, 
in a cosmology-independent way, the parameters of the 
reported correlations are given assuming an input cos- 
mology and obviously depend on the same cosmological 
parameters that we would like to constrain. Thus, apply- 
ing such relations to constrain cosmological parameters 
leads to biased results. In Ref.[47j the circular problem 
is naturally eliminated by marginalizing over the free pa- 
rameters involved in the correlations; in addition, some 
results show that these correlations do not change signif- 
icantly for a wide range of cosmological parameters [48| . 
Therefore, in this paper we use the 69 GRBs sample over 
a redshift range from z = 0.17—6.60 published in Ref.[26j 
but we keep in mind the issues related to the "circular 
problem" that are more extensively discussed in Ref. [13] . 

In Table I we list two types of linear growth factors 
data we use: i) the normalization erg inferred from the 
SDSS Lya power spectrum^ [Hjll^l; ii) the linear growth 
rate / = fl'^ from galaxy power spectrum at low redshift 
[sH and Lya growth factor measurement obtained from 



1 In this paper, we find that the low redshifts probes could not 
constrain on the early dark energy models very well. Therefore, 
we focus on the constraints on the EDE model from high red- 
shifts probes. Furthermore, in Ref. [45l| we discussed the effect of 
constraints on dark energy models with or without systematic 
error. When including the systematic error, the constraints are 
slightly weaker, but the main results are not changed. 

^ These data at different redshifts are not totally independent. The 
derived erg values are obtained by using the covariance matrix 
as measured from the data of the SDSS flux power spectrum 
and there arc (very weak) correlations in the flux power between 
different redshift bins. So, although not totally independent, we 
regard the derived values of the power spectrum amplitude from 
the SDSS flux power as nearly independent. 
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TABLE I. as values and growth rates with la error bars 
from which we derived the Unear growth factors used in our 
analysis. 
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the SDSS Lya power spectrum at z = 3 by Ref. [25| . 

We have used two Lya forest data sets: i) the high 
resolution QSO absorption spectra presented in Ref.[49[ 
consisting of the LUQAS sample [S^l and the reanalyzed 
data in Ref.jsl] (C02, this data set as a whole will be 
labelled as VHS in the following); it) the SDSS Lya for- 
est sample presented in Ref. |54| . The SDSS Lya forest 
data set consists of 3035 QSO spectra with low resolu- 
tion {R 2000) and low S/N (< 10 per pixel) span- 
ning a wide range of rcdshifts {z = 2.2 — 4.2). while 
the LUQAS and the C02 samples contain 57 high reso- 
lution (i? ~ 45000), high signal-to-noise (> 50 per pixel) 
QSO spectra with median redshifts of z = 2.125 and 
z — 2.72, respectively. The flux power spectrum of the 
Lya forest is the quantity which is observed and needs 
to be modeled at the percent or sub-percent level us- 
ing accurate numerical simulations that incorporate the 
relevant cosmological and astrophysical processes, in or- 
der to extract the underlying (linear) dark matter power 
spectrum. In this paper, we will use the derived linear 
power spectrum measured by the data set of VHS, based 
on "VHSLya" together with the SDSS Lya power spec- 
trum "SDSSLya", since they are in agreement and this 
latter has a stronger constraining power. More precisely, 
the VHSLya power spectrum consists of estimates of the 
linear dark matter power spectrum at nine values in the 
wavenumber space k at z = 2.125 and nine values at 
z ^ 2.72, in the range 0.003 < A:(s/km)< 0.03, while the 
SDSSLya consists of a single measurement at z = 3 and 
k — 0.009 s/km of amplitude, slope and curvature. The 
estimate of the uncertainty of the overall amplitude of 
the matter power spectrum is 29% for the first dataset 
and ~ 14% for the second. This estimate takes into ac- 
count possible systematic and statistical errors (see the 
relevant tables of VHS for a detailed discussion). The 
code assigns a Gaussian prior to the corresponding nui- 
sance parameter and marginalize over it. Finally, we will 
also rely on the growth factor measurements from the 



TABLE II. Assumed experimental specifications for the 
mock Planck-like measurements. The noise parameters At 
and Ap are given in units of /.tK-arcmin. 
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two data sets, derived using numerical simulations of the 
observed flux power, and we will label this as GFLya. 



2. Future Measurements 

In order to forecast future measurements we will use 
the same observables as before without BAO but exploit- 
ing the role of 21cm maps at high redshifts. 

For the simulation with PLANCK [HBl, we follow the 
method given in Ref. [5^ and mock the CMB temper- 
ature (TT) and polarization (EE) power spectra and 
temperature-polarization cross correlation (TE) by as- 
suming a given fiducial cosmological model. In Table II, 
we list the assumed experimental specifications of the fu- 
ture (mock) Planck measurement. 

The proposed satellite SNAP (Supernova / Acceler- 
ation Probe) will be a space based telescope with a 
one square degree field of view that will survey the 
whole sky [s^l- It aims at increasing the discovery rate 
of SNIa to about 2000 per year in the rcdshift range 
0.2 < z < 1.7. In this paper we simulate about 2000 
SNIa according to the forecast distribution of the SNAP 
[sst . For the error, we follow the Ref. [s^] which takes 
the magnitude dispersion to be 0.15 and the systematic 
error cTgys = 0.02 x z/1.7. The whole error for each data 
is given by CTmag(2i) = y^a^^^(z^yTol^V^ ^ where is 

the number of supernovae of the i'th rcdshift bin. Fur- 
thermore, we add as an external data set a mock set of 
400 GRBs, that mimic SWIFT ^ observations, in the 
redshift range < z < 6.4 with an intrinsic dispersion in 
the distance modulus of cr^ = 0.16 and with a redshift 
distribution very similar to that of Figure 1 of Ref. [5§] . 

For the linear growth factors data, wc simulate the 
mock data from the fiducial model of Table I with the 
error bars reduced by a factor of two. This is probably 
reasonable given the larger amounts of Lya forest data 
(e.g. SDSS-III or the X-shooter spectrograph [glj) that 
will become available soon as long with a better control 
of several systematics errors and more importantly on 
the thermal history of the intcrgalactic medium. 

We also simulate the Lya power spectrum, consisting 
of three data points at z = 2, 3 and 3.5 and at wavenum- 
bers of 0.002 s/km, 0.009 s/km and 0.02 s/km, with 1% 
fractional errors in a similar way as done by Ref. [6^ . in 
which a forecasting on the perspective of constraining the 
neutrino masses was presented. 
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The availability of 21cm (1420 MHz) maps with fu- 
ture experiments like LOFAR [gJ will allow to mea- 
sure in a given frequency interval and for a given po- 
sition in the sky the differential brightness temperature 
of a patch of gas against the CMB. Thus a completely 
new window on the neutral hydrogen content of the uni- 
verse, regarded as a tracer of cosmic structures and of 
astrophysical processes, will be opened soon. The AP 
test [231 assumes that the for an isotropically distributed 
set of astronomical objects the separations along and 
across the line-of-sights scale differently and thus gen- 
erate anisotropics. The geometrical distortion depends 
on the quantity [H{z) Da{z)]~^ that thereby can be 
measured and provide constraints on some parameters. 
The extended LOFAR will operate between 30 and 240 
MHz, with spatial resolution of the order of arcseconds, 
probing the universe in the redshift range z ^ 6 — 45. 
For our purposes here we will assume that the quantity 
[H{z) Da{z)]~^, will be measured at z = 6.5, 8, 12 with 
5% precision. Although other measurements could be 
performed using 21cm datasets, such as (quite optimisti- 
cally) a full measurement of the matter power spectrum 
at very high redshifts, we will concentrate here only on 
the AP test, since this appears to be easier to be per- 
formed and it is less sensitive to the astrophysical uncer- 
tainties that relate the differential brightness tempera- 
ture to the matter density (see however Ref. [G^I for other 
possible and very promising constraints). It is neverthe- 
less important to note that LOFAR sensitivity drops for 
z > 12 and useful HI maps at these redshifts could only 
be obtained with future ground based telescopes like SKA 
[6^ and with a better understanding of the foregrounds 
that contaminate the cosmological signal (see for exam- 
ple the discussion in Ref.[30l|). 

Here, we will use a somewhat more conservative ap- 
proach and simulate a measurement of the quantity de- 
scribed above in a redshift range accessible to the LOFAR 
experiment. Although challenging, we believe that the 
AP test could be performed and give constraints similar 
to those obtained here. 



III. RESULTS FROM PRESENT DATASETS 

In our analysis, we perform a global fitting using the 
CosmoMC package. We assume purely adiabatic ini- 
tial conditions and a flat universe, with no tensor con- 
tribution. We vary the following cosmological parame- 
ters with top-hat priors: the dark matter energy density 
rich'^ e [0.01,0.99], the baryon energy density Qbh'^ G 
[0.005,0.1], the primordial spectral index G [0.5, 1.5], 
the primordial amplitude log[10^°As] G [2.7,4.0] and the 
angular diameter of the sound horizon at last scattering 
9 G [0.5, 10]. For the pivot scale we set kso = 0.05 Mpc~^. 
When CMB data are included, we also vary the optical 
depth to reionization r G [0.01,0.8]. We do not consider 
any massive neutrino contribution. From the parame- 
ters above the MCMC code derives the reduced Hubble 
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Figure 2: ID current (marginalized) constraints on 
the dark energy parameters wo, C, as well as ag, 
wede(ziss), nEDE(2iss) and ^EDE.sf, from different current 
data combinations: WMAP5+BA0-)-SN (red solid lines), 
WMAP5+BAO+SN+GRB+Lya (blue dash-dot lines) and 
WMAP5-FBAO-hSN+GRB-HGFLyQ (black dashed lines). 



parameter Hq, the present matter fraction Jlmo, and 
f^EDE.sf: so, these parameters have non-flat priors and 
the corresponding bounds must be interpreted with some 
care. In addition, CosmoMC imposes a weak prior on the 
Hubble parameter: h G [0.4,1.0]. For part of the anal- 
ysis, we also present constraints on the growth index 7, 
defined as: 



din (5(a) 
rflna 



(6) 



This parameter, introduced by [66| . provides a good fit 
to the growth rate for many different cosmological mod- 
els and depends on the growth of perturbations, H(z) 
and f^mo- It has been recently measured by Refs.[67l.l68l| 
using mainly the observed growth factors at high and 
low redshifts and values different from that expected for 
the pure ACDM model (7 ^ 0.55) could be important 
for constraining both quintessence models and modified 
gravity scenarios (see Ref. [69l. [70j for forecasting and a 
more extensive discussion on linear growth rate in mod- 
ified gravity models). 



A. CMB+BAO+SNIa 

In Table III we show our main global constraints on 
the early dark energy model from different current data 
combinations. 

Firstly, we present the constraints from the data com- 
bination: CMB+ BAO+SNIa. In Fig. [2] we show the one 
dimensional distributions of some cosmological parame- 
ters from this data combination (solid red lines) . For the 
parameters describing the equation of state of early dark 
energy, the constraints are still weak, namely the 95% 
upper limits are wq < —0.906 and C < 2.711. The limit 
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TABLE III. Constraints on the early dark energy model from the current observations. Here we show the mean and the best 
fit values. For some parameters that are only weakly constrained we quote the 95% upper limit. 



Parameter 


WMAP5+BA0+SN 


+GRB+Lya 


+GRB+GFLya 


Mean 


BestFit 


Mean 


BestFit 


Mean 


BestFit 


Wo 


< -0.906 


-0.972 


< -0.952 


-0.999 


< -0.958 


-0.997 


C 


< 2.711 


1.858 


< 2.613 


1.628 


< 2.245 


0.152 




0.261 ± 0.014 


0.258 


0.283 ±0.014 


0.279 


0.275 ±0.014 


0.274 


0-8 


0.748 ± 0.049 


0.734 


0.842 ± 0.022 


0.863 


0.836 ± 0.024 


0.846 


^^EDe(21ss) 


< 0.0228 


0.0064 


< 0.0029 


1.77 X 10"® 


< 0.0014 


1.71 X 10"^ 


f^EDE.sf 


0.0643 ± 0.0076 


0.0672 


0.0540 ± 0.0029 


0.0529 


0.0546 ± 0.0024 


0.0543 
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0.622 ±0.139 


0.552 







on the present equation of state of dark energy compo- 
nent is still consistent with other recent results in the 
literature [J, However, current observational data, 
could not determine the "running" of equation of state 
very well. The reason is that the data we use are only 
the low redshift BAO and SNIa (z < 2), as well as the 
very high redshift CMB data at last scattering surface 
{z 1090). Due to the lack of the cosmological probes 
in the dark ages (2 < z < 1100), we do not know exactly 
the evolution of equation of state of early dark energy. 
Therefore, the constraint on the "running" of equation 
of state, which is parameterized by C, is rather poor. 

Consequently, current observations still allow very 
large amount of early dark energy at high redshift z ~ 
1090 as r2EDE(^iss) < 0.0228 (95% C.L.), which is con- 
sistent with the results obtained by Refs.fl^, [zH- Early 
dark energy models with a non-negligible fraction of dark 
energy density still fit the data very well. Furthermore, 
there are two peaks in the posterior distribution of equa- 
tion of state of early dark energy at CMB last scattering 
surface wede(2:iss)- The first peak around wede(2iss) = 
is in fact caused by the poor constraint on C. When C is 
larger than one, the equation of state wede will approach 
zero at high redshift. The larger C is, the faster wede 
approaches w — (we will come back to this point in 
Sect. IV). Considering the weak constraint on C, there 
is significant probability that allows ■u;ede(ziss) = 0. The 
second peak around wede(-ziss) = — 1 can be understood 
easily since the pure ACDM model is still favored by the 
current observations. When wq is close to —1 and C is 
close to 0, which is consistent with the data, ti'EDE(-2iss) 
reaches the cosmological constant value. Therefore, in 
order to remove the first peak, we need a much more 
stringent constraint on C from future measurements. 

As we showed in Sect. II, early dark energy models 
also affect structure formation Q. From Fig. |T] we can 
see that the linear growth factor of early dark energy 
model is suppressed significantly by the large value of 
parameter C. The non-negligible dark energy density 
at high redshifts slows down the linear growth function 
of the matter perturbation and leads to a low value of 
as today. Using CMB+BAO+SNIa data combination, 
illustrated in Fig. [21 we obtain the limit on erg today 
of as ^ 0.748 ± 0.049 (68% C.L.), which is obviously 
lower than one obtained in the pure ACDM framework: 




C «EDE<^lss) 

Figure 3: 2D (marginalized) contours of some 
parameters from different current data combi- 
nations: WMAP5±BAO±SN (red solid lines), 
WMAP5±BAO+SN+GRB+Lya (blue dash-dot lines) 
and WMAP5+BA0+SN+GRB+GFLya (black dashed 
lines) . 



as = 0.812 ± 0.026 while the error bar is enlarged 
by a factor of two. Therefore, the as today and C are 
anti-correlated, as shown in Fig. [S] On the other hand, 
the as and equation of state wq at present also have an 
anti-correlation. When wq is increased, the fraction of 
dark energy density becomes large, which also leads the 
lower as today. 

Finally, we discuss the average value of the fraction 
of dark energy in the structure formation era riEDE,sf- 
In Table III we can find that the current constraint is 
f^EDE.st = 0.0643 ± 0.0076 at la confidence level. From 
Eq.Q and Fig. [3] it is easy to sec that a large r^EDEl^iss) 
leads to the high r^EDE.sf- In Fig. [4] we also plot the 
two dimensional plots between f^EDE.sf and other related 
parameters. As we discussed before, when wq and C 
become large, the fraction of dark energy density at high 
redshifts will also be large. On the other hand, when 
f^EDE.sf increases, the linear growth function of matter 
perturbation will be slowed down and the value of as 
today will decrease consequently: riEDE,sf and the as 
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Figure 4: 2D (marginalized) contours of S^EDE.sf ver- 
sus wo, C, (78 and iyEDE(2:iBB) from different current 
data combinations: WMAP5-I-BAO+SN (red solid lines), 
WMAP5+BAO+SN+GRB+Lya (blue dash-dot lines) and 
WMAP5+BA0+SN+GRB+GFLya (black dashed lines). 



today arc then obviously anti-corrclated. 

To sum up, using CMB, BAO and SNIa, we constrain 
the amount of dark energy at the last scattering surface 
to be less than 2% and in the structure formation era 
to be of the order of 6%. This last number is in rough 
agreement with the results of [zl], who used WMAP 3 
year data and some external data sets, although their 
values are somewhat smaller (~ 4%) and this is due to 
the slightly different definition of SlEDE.sf • 



B. Adding High-Redshift Probes 

From the results above, we can see that without ob- 
servations in or around the so called dark ages, the con- 
straints on the early dark energy model are interesting 
but not very tight. In this subsection we add some 
high-redshift observational data, such as GRBs, Lya and 
GFLya data, to improve our constraints. Because most 
of GFLya data are obtained from the Lya forest power 
spectrum, in order to avoid overestimating the weight 
of Lya data, in our analysis we do not use GFLya and 
Lya data together. And the results, as we expect, from 
these two data combinations arc consistent between each 
other. Therefore, in the following we mainly present the 
global fitting results obtained from those that include 
GRB and GFLya data. 

When adding the high-redshift probes, the constraints 
on and C improve significantly: wo < —0.958 and 
C < 2.245 at 95% confidence level, while their best fit 
values are very close to the pure ACDM model. Natu- 
rally, the 95% upper Hmit of I^edeI^Iss) has also been 
tightened by a factor of ten, riEDE(^iss) < 1-4 x 10""^. A 
large number of early dark energy models have been ruled 



out by the high-redshift probes. The high-redshift obser- 
vations are indeed effective in constraining further the 
early dark energy models investigated here. Note that 
adding only the power spectrum information as provided 
by Lya improves the constraints on fiEDE(2iss) by a fac- 
tor five, when compared to CMB-|-BAO-|-SNIa, so the 
growth factor information is really fundamental in con- 
straining the EDE model. However, the double peak in 
the distribution of wedeI-^Iss) is still present, which im- 
plies that the current data are not accurate enough to 
measure wede('ZIss), even when high-redshift probes are 
added. 

Another effect of adding high-redshift probes is the 
larger value of erg today. As we know, the Lya forest 
data favor a larger erg ^ 0.9 than the CMB measure- 
ment'^. When we add the GFLya or Lya forest data, we 
obtain a higher erg today than before: erg ~ 0.836 ± 0.024 
(Icr). However, this is still smaller than the value favored 
by Lya forest data alone in the pure ACDM model and 
consistent with WMAP5-I-BA0-I-SN measurement of 
We also note that if we use only VHSLya data sets and 
adding them to WMAP5+BA0+SN and GRBs the re- 
sults do not change significantly as compared to having 
SDSSLya and VHSLya together. This is mainly due to 
the constraining power of the growth factors measure- 
ment obtained with VHSLya. 

The conclusion of this section is that adding higher red- 
shift probes constrains the energy density of dark energy 
at the last scattering surface to be around 0.1%. 



C. Growth Index 7 

Finally, we extend our discussion by briefly addressing 
the perspective of constraining the linear growth index 7. 
Measurements of the growth history of cosmic structures 
combine information on both cosmic expansion and the 
underlying theory of gravity, the parameter 7 is thus a 
unique prediction of any modified gravity model and a 
simultaneous fitting of this parameter and other dynam- 
ical probes could provide constraints on modified gravity 
scenarios and dynamical dark energy models. 

By using all the current data available we obtain 
7 = 0.622 ± 0.139 (Ict) error bar for our EDE model, 
which is in agreement with the ACDM and with the 
values obtained by Refs.[63, [6^ (with slightly smaller 
error bars). The present error on 7 is also similar to 
that inferred from future weak lensing and SNIa data by 
Refs.m. 



^ This fact docs not imply a contradiction or a failure of the ACDM 
model but there is indeed a tension between small scale con- 
straints and the large scale ones. However, more recent stud- 
ies that combine WMAP year 3, the Lyman-a forest and weak- 
lensing measurements of the COSMOS-z survey [2^ point to a 
value of (Tg = 0.800 it 0.023, which is in good agreement with the 
findings of WMAP year 5 plus BAO and SNe 
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Figure 5: The posterior distribution of 7 and 2D (marginal- 
ized) contours of 7 and wo, C, wede(2iss), nEDE(ziss) and 
^EDE.sf from WMAP5+BAO+SN+GRB+Lya (blue dotted 
lines) . 
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Figure 6: ID current (marginalized) constraints on the dark 
energy parameters wo, C, as weU as as, u;ede(ziss), nEDE(aiss) 
and f^EDE.sf from future measurements with the fiducial mod- 
els; ACDM (red dashed fines), EDEl model (black sofid 
lines) . 



TABLE IV. Constraints on the early dark energy model 
from future measurements. We also list the standard 
deviation of these parameters based on the future mock data 
sets. For parameters that are only weakly constrained we 
quote the 95% upper limit instead. 



Parameter 


Fiducial ACDM 


Fiducial EDE 


Wo 


< -0.975 


[-0.985, -0.920] 


C 


< 2.635 


[0.256, 2.300] 




0.0029 


0.0032 


as 


0.0078 


0.0145 


f2EDE(2lss) 


< 1.02 X 10"'* 


< 0.011 


f^EDE.sf 


0.0011 


0.0037 


7 


0.026 


0.048 



In Figurc[5]we show the 2-dimcnsional contours of 7 vs. 
CT8, C*, Wo and ilEDE(^iss) to address possible degeneracies 
of this parameter. For the current data sets, we note that 
the degeneracies are weak and only the degeneracy in the 
7 — f2EDE,sf plane is significant. This is easily understood 
since a larger 7 implies a faster growth of structures and 
thereby a smaller amount of dark energy. 



IV. RESULTS FROM FUTURE DATA SETS 

Since the present data do not give very stringent con- 
straints on the parameters of early dark energy model, 
it is worthwhile discussing whether future data could de- 
termine these parameters conclusively. For this purpose 
we have performed a further analysis and we have cho- 
sen two fiducial models in perfect agreement with current 
data: a pure ACDM model and an EDE model with pa- 
rameters taken to be the best-fit values of Table III from 
the current constraints of CMB+BAO+SNIa. This latter 
model will be labelled EDEl, and its evolution is shown 
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Figure 7: Evolution of different equation of state for the 
Mocker model when fixing C = 2. All these models fit low 
redshift (z < 2) and CMB constraints but have very different 
evolution at intermediate redshifts. 



in Figdl 



A. Fiducial ACDM 

Firstly, wc choose the pure ACDM as the fiducial 
model. In Table IV we list the forecasts for some pa- 
rameters using all the future measurements described in 
Section |TTB2j GRBs, SNIa, Lya (matter power spec- 
trum and growth factors) and the AP test. 

Due to the smaller error bars of the mock data sets, 
the constraint on wq improves significantly, namely the 
95% upper limit is now wq < —0.975. The fraction of 
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Figure 8: 2D marginalized contours C vs Wit from future mock 
measurements of the fiducial EDEl model. 



early dark energy density has been limited very stringent, 
f^EDE(ziss) < 1.02 X 10-'' (95% C.L.), which is reduced 
by another factor of ten, compared with the present con- 
straint. The error on the 7 parameter that we forecast 
is 0.026, while the error on the amount of EDE in the 
structure formation era is about 1.1 x 10~^. Therefore, 
we can foresee a very precise determination of the amount 
of dark energy in the dark ages (or in the structure for- 
mation era) with this future data set. 

By contrast, the 95% upper limit on C has not been im- 
proved, C < 2.635. This is due to the chosen parameteri- 
zation of early dark energy. In Fig. [7] we illustrate the dif- 
ferent evolution of several EDE equations of state varying 
wq and keeping C = 2. We note that the smaller the red- 
shift at which wq approaches the cosmological constant, 
the smaller the redshift at which w{z) flattens to zero in 
the high redshift universe. For example, if wq = —0.95, 
we have w{z = 10) ~ —0.1 and C = 2 has been ruled 
out obviously. On the other hand, if wq = —0.9999, we 
still have w{z = 10) ~ —1. Consequently, C = 2 is still 
favored by the data. Therefore, even if the constraint on 
Wq becomes stringent, a large value of C is still allowed, 
and then the constraint on C will be somewhat poorer. 
Moreover, due to the allowed large value of C, the first 
peak of wede(-Ziss) around = is still present, as shown 
in Fig. [51 where the ID marginalized constraints at the 
last scattering surface are shown. 



B. Fiducial EDE 

For our EDE fiducial model we choose the best fit val- 
ues of current constraint as the fiducial model to simulate 
the future measurements, wq = —0.972 and C = 1.858 
(EDEl). 

In Fig. [SI we plot the one dimensional distributions 
of Wq and C. The 95% confidence level is —0.985 < 
Wo < -0.920 and 0.256 <C < 2.300, respectively. These 



results imply that the future measurements could dis- 
tinguish the pure ACDM model and early dark energy 
models at least at 5a confidence level. We also show 
the contour in the [wojC] plane in Fig. [51 We note that 
even with these futuristic data sets simulated in a con- 
servative way the value C = is still allowed at a 2a 
confidence level. When wq increases, C must decrease in 
order to match the current observations. Thus, wq and 
C are anti-corrclated. Consequently, the allowed value of 
f^EDE(2iss) could be very large: f^EDE(^iss) < 1-1 x 10"^ 
(95% C.L.). Furthermore, with this fiducial model, the 
double peak in the distribution of ?i'EDE(2iss) disappears: 
because the pure ACDM is disfavored at more than 5a 
confidence level. 

In Fig.[6lwe clearly appreciate the very distinctive pre- 
dictions of ilEDE.sf and fJs for the ACDM and EDEl mod- 
els. The values obtained for l^EDE.sf and ag are different 
by 20% and 10%, respectively, and can be measured at 
very high precision in the two cases. 

We are also interested in understanding which of the 
different cosmological probes give the most stringent con- 
straints. In order to investigate this we perform three 
different runs: we find that the best constraints are ob- 
tained when future Lya and GFLya are added to Planck 
and SNAP data, while adding AP only to Planck and 
SNAP results in the the worst constraints. In the third 
run, we consider adding only GRBs mock data to Planck 
and SNAP and in this case the results are in between 
those inferred from Lya and AP test. The degeneracies 
between different parameters are similar in three cases. 
We also checked that degrading the accuracy of the AP 
test to 10% as opposed to the chosen 5% has not a sig- 
nificant impact on the recovered parameters. 

We also consider the results in terms of 7 has done in 
Section IIII CI We find that we could constrain 7 to high 
significance in the two cases: the standard deviations are 
0.026 and 0.048, which are reduced by a factor of five 
and three, when compared to the current constraints, 
respectively. By looking at Fig. 9, we can now see that 
the constraints on 7 are more precise and degeneracies 
with Oede(2;iss) and r^EDE.sf are stronger. 



V. CONCLUSIONS AND DISCUSSION 

With the availability of new data from SNIa, and the 
advent of relatively new cosmological probes, such as 
Lya and GRBs in the high redshift universe, there has 
been an increasing interest in the study of dynamical dark 
energy models with a significant energy contribution in 
the structure formation era. In this paper we present con- 
straints on a particular early dark energy model using the 
latest observations. We use the Mocker parameterization 
proposed by Ref . [Ill , but our results are general since the 
amount of dark energy at the last scattering surface and 
in the structure formation era have also been quoted in 
order to better cornpare with other possible dark energy 
parameterizations [72l |. 
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Figure 9: The posterior distribution of 7 and 2D (marginal- 
ized) contours of 7 and wo, C, wede(2iss), nEDE(zias) and 
nEDE,sf from future measurements with the fiducial models: 
ACDM (red dashed lines), EDEl model (black solid lines). 



We find that current constraints based on CMB, BAO 
and SNIa data limit the amount of dark energy density at 
the last scattering surface to be 2%, while in the structure 
formation period is of the order of 6%. If we add high 
redshift GRBs and the Lya forest data, in the form of a 
matter power spectrum or in the form of a growth factor 
measurement, the constraints improve by a factor five 
and fifteen, respectively. The ACDM model is still a very 
good fit to the data currently available and only upper 
limits on the amount of dark energy at the last scattering 
could be quoted. 

We also forecast what could be done with future data 
sets in particular: Planck for the CMB; GRBs and SNIa 
luminosity distances as could be measured by SNAP; a 
better measurement of Lya matter power spectrum and 
growth factors as could be obtained from SDSS-III and 
the AP test performed on HI 21 cm maps from LOFAR. 
With these futuristic mock data sets we simulate two 
different fiducial models that are currently in agreement 
with observations a ACDM model and an EDE model 
with about 6% contribution to the energy density in the 
structure formation epoch. We find that the two mod- 
els will be clearly distinguishable at more than 5a con- 
fidence level. Thereby, the role of these high redshift 
probes is crucial in discriminating between pure ACDM 
and quintessence and/or modified gravity models. 

Furthermore, motivated by the fact that the growth 
factors at high redshift seem to have the highest con- 
straining power, we also address the perspectives of mea- 
suring the parameter 7 that parameterizes the growth 
rate of structures and is usually a very precise predic- 



tion for any quintessence and/or modified gravity model. 
Even in this case, we find that in the two models we could 
constrain 7 to high significance: the standard deviations 
are reduced by a factor of five and three, when compared 
with the current constraint, respectively. 

We conclude this Section with some words of cau- 
tion: while high redshift probes are important in con- 
straining early dark energy models, the obtained values 
from present and future constraints rely mainly on the 
Lya forest and its derived growth factors. These mea- 
surements, summarized partly in Table I, have been ob- 
tained in the framework of pure ACDM model and are 
valid only when small deviations from this model are con- 
sidered. We believe that we are in such a situation here, 
where the amount of dark energy has little impact on 
observed Lya quantities and this effect is possibly de- 
generate with other parameters that have been already 
marginalized over to obtain the final estimate. However, 
further work using both numerical hydrodynamical sim- 
ulations as long with an accurate control of systematics 
should be done in order to better check this (see for ex- 
ample Ref.|73j]). Along the same lines, while the role 
of SNIa as cosmological probes is somewhat reasonably 
well established, those of GRBs is still controversial. We 
tried to circumvent this problem by marginalizing over 
the nuisance parameters that enter the calibration of the 
GRBs distance moduli but, even in this case, more obser- 
vational and theoretical work is needed in order to trust 
quantitatively the results. Even regarding the forecast- 
ing using the AP test, we tried to implement these con- 
straints in a conservative way, in a range of redshift that 
could be observed by LOFAR, without assuming that a 
full matter power spectrum could be measured from the 
observed brightness temperature (this of course will have 
a fundamental impact especially at even higher redshifts 
than those considered here). 

However, even considering these cautious remarks, our 
work clearly emphasizes the capabilities of present and 
future data in an intermediate redshift range which is 
poorly probed by the observations, to quantitatively con- 
strain the dark energy component and its evolution. We 
consider these observations as very important since they 
can open up a completely new window on observational 
cosmology and high energy physics and provide new in- 
sights on cither modified gravity or quintessence scenar- 
ios. 
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